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Abstract. Stability of Wardrop equilibria is analyzed for dynamical transportation networks 
in which the drivers' route choices are influenced by information at multiple temporal and spatial 
scales. The considered model involves a continuum of indistinguishable drivers commuting between 
a common origin/destination pair in an acyclic transportation network. The drivers' route choices 
are affected by their, relatively infrequent, perturbed best responses to global information about the 
current network congestion levels, as well as their instantaneous local observation of the immediate 
surroundings as they transit through the network. A novel model is proposed for the drivers' route 
choice behavior, exhibiting local consistency with their preference toward globally less congested 
paths as well as myopic decisions in favor of locally less congested paths. The simultaneous evolution 
of the traffic congestion on the network and of the aggregate path preference is modeled by a system 
of coupled ordinary differential equations. The main result shows that, if the frequency of updates 
of path preferences is sufficiently small as compared to the frequency of the traffic flow dynamics, 
then the state of the transportation network ultimately approaches a neighborhood of the Wardrop 
equilibrium. The presented results may be read as a further evidence in support of Wardrop's 
postulate of equilibrium, showing robustness of it with respect to non-persistent perturbations. The 
proposed analysis combines techniques from singular perturbation theory, evolutionary game theory, 
and cooperative dynamical systems. 
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1. Introduction. As transportation demand is dramatically approaching its in- 
i— i frastructure capacity, a rigorous understanding of the relationship between the macro- 

scopic properties of transportation networks and realistic driver route choice behavior 
is attracting renewed research interest. Such an analysis is essential, among other 
things, for appropriate design of incentives influencing drivers' behavior in order to 
(N) induce a desired socially optimal usage of the transportation infrastructure. A partic- 

ularly relevant issue is the impact of drivers' en route responses to unexpected events 
on the overall transportation network dynamics. This issue is particularly signifi- 
cant in a modern real-life transportation network scenario, where recent technological 
advancements in intelligent traveller information devices have enabled drivers to be 
much more flexible in selecting their routes to destination even while being en route. 
While there has been a significant research effort to investigate the effect of such tech- 
nologies on the route choice behavior of drivers, e.g., see [19, 16], the analytical study 
of the dynamical properties of the whole network under such behavior has attracted 
5_i very little attention. 
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This paper is focused on the stability analysis of transportation networks in a 
setup where the drivers have access to traffic information at multiple temporal and 
spatial scales and they have the flexibility to switch their route to destination at 
every intermediate traffic intersection. Specifically, we consider a model in which 
the drivers choose their routes while having access to relatively infrequent global 
information about the network congestion state, and real-time local information as 
they transit through the network. The drivers' route choice behavior is then influenced 
by relatively slowly evolving path preferences as well as myopic responses to the 
instantaneous observation of the local congestion levels at the intersections. This 
setup captures many real-life scenarios where unexpected events observed en route 
might cause drivers to take a temporary detour, but not necessarily to change their 
path preferences. Such path preferences may instead be updated, e.g., on a daily, 
weekly, or longer time basis, in response to information about the global congestion 
state of the different origin-destination paths collected from the drivers' personal 
experience, their opinion exchanges with their peers, as well as from information 
media. However, since the traffic dynamics is significantly influenced by the drivers' 
response to real-time local information, such responses can influence the drivers' path 
preference thereby modifying their global route choice behavior in the long run. We 
propose and analyze a novel model for the drivers' route choice behavior that combines 
relatively infrequent information about the global congestion status of the network 
with real-time local observations as explained below. 

In our model, the network is represented by a directed acyclic graph with one 
origin and one destination. A continuous constant flow of indistinguishable drivers 
enters from the origin, and flows through the network until reaching the destination 
node. Traffic parameters, such as average speed, traffic density, and flow, are modeled 
as homogeneous quantities on every link, related one to each other by functional 
dependencies representative of the links' congestion properties. The dynamics of 
such traffic parameters is governed by the law of conservation of mass, as well as 
the drivers' route choice behavior. In turn, the drivers' route choice behavior is 
assumed to be influenced by two factors: the aggregate path preference, measuring 
the relative appeal of the different routes to the drivers, and local observations of the 
current congestion levels. The path preference dynamics evolve at a slow time scale 
(as compared to the traffic dynamics), following a perturbed best response to global 
information, embodied by the current congestion levels on the whole network. When 
traversing an intermediate node in the network, drivers behave according to their path 
preference, if this is consistent with the current, locally observed, aggregate behavior 
of the other drivers. On the other hand, when there is a discrepancy between the 
aggregate path preference and the locally observed aggregate behavior, then drivers 
tend to compensate this by myopically preferring routes which appear to be locally 
less congested. 

The above-described model gives rise to a double feedback dynamics, governed by 
a finite-dimensional system of coupled ordinary differential equations. Such a dynam- 
ical system has two natural time scales, characterizing the dynamics of the drivers' 
aggregate path preference and of the traffic parameters on the different links, respec- 
tively. We study the long-time behavior of this dynamical system: our main result 
shows that, in the limit of small update rate of the aggregate path preferences, a state 
of approximate Wardrop equilibrium [20] is approached. The latter is a configuration 
in which the delay associated to any source-destination path chosen by a nonzero 
fraction of the drivers does not exceed the delay associated to any other path. Our 
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results contribute to providing a stronger evidence in support of the significance of 
Wardrop's postulate of equilibrium for a transportation network. They may also be 
read as a sort of robustness of such equilibrium notion with respect to non-persistent 
perturbations of the network. 

The analytical arguments wc propose mainly rely on three ideas: adopting a 
singular perturbation approach [13], by considering the aggregate path preference as 
'quasi-static' when studying the fast scale dynamics of the traffic parameters, and the 
traffic parameters as 'almost equilibrated' when analyzing the slow scale dynamics of 
the aggregate path preference; exploiting the inherent cooperative 1 dependence of the 
route choice function on the local traffic parameters in order to establish exponential 
stability of the fast scale dynamics of the traffic parameters; adapting results from 
evolutionary population games [12, 18] in order to establish stability properties of the 
slow scale perturbed best response dynamics of the aggregate path preference. 

Our work is naturally related to two streams of literature on transportation net- 
works. On the one hand, traffic flows on networks have been widely analyzed with 
fluid-dynamical and kinetic models: see, e.g., [8], and references therein. As compared 
to these models (typically described by partial, or intcgro-diffcrential equations), ours 
significantly simplifies the evolution of the traffic parameters (treating them as homo- 
geneous quantities on the links, representative of spatial averages), whereas it high- 
lights the role of the drivers' route choice behavior with its double feedback dynamics, 
which is typically neglected in that literature. 

On the other hand, transportation networks have been studied from a decision- 
theoretic perspective within the framework of congestion games [3, 17]. In these 
models, drivers make sequential myopic route choice decisions in pursuit of mini- 
mizing their personal travel times, in response to complete information about the 
whole network. Congestion games are known to belong to a class of games known 
as potential games, a consequence of which is that, best responses of the drivers are 
aligned with the gradient of a common potential function and hence the system even- 
tually converges to a critical point of this potential function, which, under appropriate 
monotonicity conditions of the congestion properties of the links of the network, cor- 
responds to a Wardrop equilibrium. Such an approach has been used, for example 
in [14]. The stability of Wardrop equilibrium in the context of communication net- 
works has been studied in [4]. It is important to note that the two salient features 
of a typical congestion game setup are that information is available to the drivers at 
a single temporal and spatial scale, and that the dynamics of traffic parameters are 
completely neglected by assuming that they are instantaneously equilibrated. In con- 
trast, we study the stability of Wardrop equilibrium in a setting where the dynamics 
of the traffic parameters are not neglected, and the drivers' route choice decisions are 
affected by, relatively infrequent global information, as well as their real-time local 
information as they transit through the network. As a consequence, classic results of 
evolutionary game theory and population dynamics [12, 18] are not directly applicable 
to our framework, and novel analytical tools have to be developed, particularly for 
the analysis of the fast scale dynamics of the traffic parameters. 

The rest of the paper is organized as follows. In Section 2, we formulate the model 
and state the main result. Section 3 is a technical section that contains the proofs for 
the main result including intermediate results. In Section 4, we report results from 
illustrative numerical experiments. Finally, we conclude in Section 5 and also mention 
potential future research directions. 



Here, the adjective 'cooperative' is intended in the sense of Hirsch [9, 10]. 
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Before proceeding, we establish here some notation to be used throughout the 
paper. Let R be the set of reals, R + :— {x € R : x > 0} be the set of nonnegative 
reals. Let A and B be finite sets. Then, \A\ will denote the cardinality of A, R A 
(respectively, R+) the space of real- valued (nonnegative-real- valued) vectors whose 
components are labeled by elements of A, and R- 4xB the space of matrices whose real 
entries labeled by pairs of elements in Ax B. The transpose of a matrix M € R AxB , 
will be denoted by M' € R Bx - 4 , while / to be an identity matrix, and 1 the al- 
one vector, whose size will be clear from the context. The simplex of probability 
vectors over A will be denoted by S(A) := {x e R^ : J2aeA Xa = 

1}. If B C A, 

1b ■ A — > {0, 1} will stand for the indicator function of B, with 1b{o) = 1 if a £ B, 
1b (a) = if a e A\B. Forp e [l,oo], || • || p is thep-norm. By default, let || -|| := || • || 2 
denote the Euclidean norm. Let int(A') be the interior of a set X C R d , and dX denote 
its boundary. Let sgn : R — > { — 1,0,1} be the sign function, defined by sgn(:r) is 1 
if x > 0, sgn(a;) = —1 if x < 1, and sgn(ir) = if x = 0. Conventionally, we shall 
assume the identity d|x|/dx = sgn(:r) to be valid for every x € R, including x = 0. 

2. Model formulation and main result. In this section, we formulate the 
problem and state the main result. In our formulation, we represent the dynamics of 
the traffic and the route choice behavior on a transportation network as a system of 
coupled ordinary differential equations with two time scales representative of route 
choice behavior influenced by the two levels of information. The key components of 
our model are: network topology, congestion properties of the links, path preference 
dynamics, and node-wise route choice decision. We next describe these components 
in detail. 

2.1. Network characteristics. Let the topology of the transportation network 
be described by a directed graph (shortly, di-graph) Q = (V, £), where V is a finite 
set of nodes and £ C V x V is the set of (directed) links. For every node v € V, 
we shall denote by £~ , and £+ , the sets of its incoming, and, respectively, outgoing 
links. A lcngth-Z (directed) path from u € V to v € V is an i-tuple of consecutive links 
{(vj_i, Vj) € £ : 1 < j < 1} with u = u, and v\ = v. A cycle is path of length I > 1 
from a node v to itself. Throughout this paper, we shall assume that: 

Assumption 1. The di-graph Q contains no cycles, has a unique origin (i.e., 
some v e V such that £~ = and a unique destination (i.e., v e V such that 
£y = $). Moreover, there exists a path to the destination node from every other node 
in V. 

Assumption 1 implies that one can find a (not necessarily unique) topological ordering 
of the node set V (see, e.g., [6]). We shall assume to have fixed one such ordering, 
identifying V with the integer set {0,1,..., n}, where n := |V| — 1, in such a way that 

C cM < £^ V« = 0,...,n. 

We shall model the traffic parameters as time-varying quantities which are ho- 
mogeneous over each link of the network. Specifically, for every link e G £, and time 
instant t > 0, we shall denote the current traffic density, and flow, by p e (t), and f e {t), 
respectively, while 

p(t) := {p e (t) : e E £} , f(t) := {f e (t) : e e £} 

will stand for the vectors of all traffic densities, and flows, respectively. Current traffic 
flow and density on each link are related by a functional dependence 



fe = Ve(Pe), ee£. 
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(2.1) 



Such functional dependence models the drivers' speed and lane adjustment behavior 
in response to traffic density on a particular segment of a road, ft will be assumed to 
satisfy the following: 

ASSUMPTION 2. For every link e e £, the flow- density function fi e : R + -> R + 
is continuously differentiable, strictly increasing, strictly concave and is such that 
fj, e (0) = and lim^o d^ e /dp e (p e ) < +00. 

Remark 1. Flow- density functions commonly used in transportation theory typ- 
ically are not globally increasing, but rather have a H-shaped graph [8j: ji e {p e ) in- 
creases from /i e (0) = until achieving a maximum C e = p, e (p e ), and then decreases 
for p e > p e . Assumption 2 remains a good approximation of this setting, provided 
that p e stays in the interval [0,p e ). 
For every link e € £, let 

C e := sup{p e (p e ) : p e > 0} = lim p e (p e ) 

be its maximum flow capacity. Moreover, let 

~ x ee£ +[0,C e ) , T := x ee£ [0,C e ) 

be the sets of local, and, respectively, global admissible flow vectors. Observe that 
our formulation allows for both the cases of bounded and unbounded maximum flow 
capacities. As the flow f e is the product of speed and density, it is natural to introduce 
the delay function 



^[0,+oof, T e (f e ):=l 



+OO if / e > C e 

p-Hfe)/fe if / e e(0,C e ), (2.2) 
I Vfe(0) if fe = 0, 



whose components measure the flow-dependent time taken to traverse the different 
links. 2 

Example 1. A flow- density function that satisfies Assumption 2 is given by 

Me(Pe) = C e (l-e- e ^) Vee£, (2.3) 
where C e > 0, and 9 e > 0. The corresponding delay function is 

T e (f e ) = ^ T log 7 ^ T . 

Veje W — Je 



We shall denote by V the set of distinct paths in Q from the origin node to the 
destination node n, and let 

be the link-path incidence matrix of Q. The relative appeal of the different paths to 
the drivers will be modeled by a time- varying probability vector over V, which will 
be referred to as the current aggregate path preference, and denoted by n(t). If one 



2 Here it has implicitly been assumed, without any loss of generality, that all the links are of unit 
length. 
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assumes, as we shall do throughout this paper, a constant unit incoming flow in the 
origin node, it is natural to consider the vector 

r ■■= A* 

of the flows associated to the current aggregate path preference. Indeed, /J = 
J2 p A ep TT p represents the total traffic flow that a link e G £ would sustain in a hypo- 
thetic equilibrium condition in which the fraction of drivers choosing any path peP 
is given by ir p . Now, let 

II := {tt G S{V) : (A7r) e < C e , Ve e £} 

be the set of feasible path preferences. Here, the term 'feasible' refers to the fact that 
the flow vector associated to any ir G II satisfies the capacity constraint /J < C e 
for every e G £. Observe that, whenever C e > 1 for every e G £ (or, in particular, 
when link capacities are infinite), the set of admissible path preferences II coincides 
with the whole simplex S{V). In contrast, when C e < 1 for some e G £, II C S(V) is 
a strict inclusion. On the other hand, whether II is empty or not depends solely on 
the value of the min-cut capacity of the network [1, Ch. 4] 

C*:= mm C u , C u := £ C e , 

OelT,niu e=(u,v)e£ : 

ueu,vev\u 

as shown in the following, simply established, result. 

Proposition 2.1. The set II is nonempty if and only if C* > 1. 

Proof. Fix a cut-set U QV such that eW, and n <£U. Then, every path peP 
contains exactly one link (u, v) € p such that u gW, and d€ V\W. Hence, for every 
7r G n, one has that 

= ^2 c e > ^epTp = ^ P = 1 ■ 

e=(«,t))s£: pGP e=(u,i))6£: peP 

ueu.vev\u ueu,vev\u 

Minimizing over all cut-sets U shows that C* > 1 is necessary for n to be nonempty. 

For the inverse implication, consider a network with the same topology Q and link 
capacities c e — max{C e — |£| _1 (C* — 1),0}. The min-cut capacity of this network 
satisfies c* > C* — (C* — 1) = 1. Note also that, from our construction, C e > c e > 0. 
Therefore, the max- flow min-cut theorem (see, e.g., [1, Thm. 4.1]) implies that there 
exists some jt £ II, thus proving that n is nonempty. ■ 

In the case when C* < 1 it is not hard to show that, since the incoming flow 
exceeds the outgoing flow of the network, the system will grow unstable, i.e., p e {t) is 
unbounded as t grows large, for some link e G £. Therefore, throughout this paper 
we shall confine ourselves to transportation networks satisfying: 

Assumption 3. The min-cut capacity satisfies C* > 1. 

2.2. Route choice behavior and traffic dynamics. We now describe the 
drivers' route choice behavior and traffic dynamics on the network. We envision a 
continuum of indistinguishable drivers traveling through the network. Drivers enter 
the network from the origin node at a constant unit rate, travel through it, and 
leave the network from the destination node n. While inside the network, drivers 
occupy some link e G £. The time required by the drivers to traverse link e, and 
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the current flow on such link are governed by its congestion properties, as given by 
(2.2), and (2.1), respectively. When entering the network from the origin node v = 0, 
as well as when reaching the tail node v € {1, 2, . . . ,n — 1} of some link e ^ £~ , 
the drivers instantaneously join some link e € £„. In this paper, we shall model the 
choice of such new link to depend on infrequently updated perturbed best responses 
of the drivers to global information about the congestion status of the whole network 
as well as on their instantaneous observation of the local congestion levels. We next 
describe these two aspects of the model in detail. 

Aggregate path preference dynamics. The drivers' aggregate path prefer- 
ence 7r(t), already introduced in Sect. 2.1, models the relative appeal of the different 
paths to the drivers' population. The aggregate path preference ir(t) is updated as 
drivers access global information about the current congestion status of the whole 
network. This occurs at some rate n > 0, which will be assumed small with respect 
to the time-scale of the network flow dynamics. Information about the current sta- 
tus of the network is embodied by the current traffic flow vector f(t). From f(t), 
drivers can evaluate the vector A'T(f(t)), whose p-th component X) e e£ A ep T e (f e (t)) 
coincides with the total delay one expects to incur on path p assuming that the con- 
gestion levels on such path won't change. Drivers' are assumed to react to such global 
information by a perturbed best response 

F h (f) := argmin {uj'A'T(f) + h(uj)} , (2.4) 

where h : Tlh — > M is an admissible perturbation, satisfying the following: 

Assumption 4. An admissible perturbation is a function h : — > R where 

Uh C II is a closed convex set, h(-) is strictly convex, twice differ -entiable in mt(Tlh), 

and is such that lim^^gUh \ \Vh(n)\\ = +00, where V := (/ — |7 , | _1 11')V is the 

projected gradient on S(P) 3 . 

As a result, the aggregate path preference ir(t) evolves as 

±* = r,{F h {f)-w) . (2.5) 

The perturbed best response function F h (f) provides an idealized description of the 
behavior of drivers whose decisions are based on inexact information about the state 
of the network. In particular, it can be shown that the form of F h (f) given in (2.4) is 
equivalent to the minimization over paths p € V of the expected delay X^ee£ A ep T e (f e ) 
corrupted by some (admissible) stochastic perturbation (see e.g. [11]). 

It is easy to establish that the perturbed best response F h (f) is continuously 
diffcrentiable on T . Moreover, it is well known [18] that, as ||/i||oo I 0, and f II, 
the perturbed best response F h (f) converges to the set argmin{u/A'T(/) : w € 11} 
of best responses. 4 

Example 2. Assume that C e > 1 for all e € £. Then, an example of perturbed 
best response satisfying Assumption 4 is the logit function with noise level /3 > 0, 
which is defined as 

Fp(/) -E 9e .exp(^T(/)),)' Per - (2 ' 6) 



3 We shall use the notation $ := / — \V\ 1 11' £ K Vx ' p to denote the corresponding projection 
matrix. 

Here, the convergences n h f n, and {F h (f)} -> argmin{VA'T(/) : u 6 It} are intended to 
hold in the HausdorfT metric, (see, e.g., [2, Def. 4.4.11]) 



This corresponds to the admissible perturbation function h(co) = —f3 1 X^ w pl°g w p- 
For any fixed f £ T , one has that linig_ ) ._|_ (x) F h {f), with F h (f) as defined in (2.6), 
is a uniform distribution over the set argmin{(AT(/)) p : p £ V}. We refer the 
reader to [7] for more on the connection between F h characterized by Assumption 4 
and smooth best response functions. 

Remark 2. In the evolutionary game theory literature, e.g., see [12, 18], the do- 
main of an admissible perturbation function h, as well as the one of the minimization 
in the right-hand side of (2.4), is typically assumed to be the whole simplex S{V), 
instead of a closed polytope 11^ C II C S(V). Notice that, as already observed in 
Sect. 2.1, when C e > 1 for every e £ £, II = S(V) is a closed polytope, so that one 
can choose Tlh = II. Therefore, in this case, Assumption 4 does not introduce any 
additional restriction with respect to such theory. 

On the other hand, when C e < 1 for some e £ £ , then the inclusions of H h C 
II C S(V) are both strict, so that Assumption 4 does introduce additional restrictions 
on the admissible perturbations. However, it is worth observing that, in a classic 
evolutionary game theoretic framework, the dynamics of the aggregate path preference 
would be autonomous rather than coupled to the one of the actual flow. In particular, 
perturbed best response dynamics in that framework would read as 

±n = F h (n-n, (2.7) 

rather than as in (2.5). For such dynamics, the fact that T e ((Air) e ) = +oo whenever 
(Air) e > C e , can be shown to imply that Tr(t) would reach a compact C II in some 
finite time and never leave it. In contrast, in the two time-scale model of coupled 
dynamics considered in this paper (see (2.11)), such more restrictive assumption is 
needed in order to ensure the same property for the trajectories of ir(t) (see Lemma 
3.4). 

Local route decisions. We now describe the local route decisions, characteriz- 
ing the fraction of drivers choosing each link e £ £+ when traversing a non-destination 
node v. Such a fraction will be assumed to be a continuously differentiable function 
Gg(/ £ +,7r) of the local traffic flow f £ + :— {f e : e £ £+}, as well as of the current 
aggregate path preference tt. We shall refer to 

G":f,xII-> S(£+) (2.8) 

as the local decision function at node v £ {0, 1, . . . , n— 1}, and assume that it satisfies 
the following: 

Assumption 5. For all < v < n, and tt £ II, 

(E je£ + ft) G ^ (fh < ^ ^ - Ve e ^ ■ 



Assumption 6. For all < v < n, tt e IT, and f £ + £ T v , 

d_ 

We 



G*(/ £ +,7r)>0, Vj/eeC. 



Assumption 5 is a consistency assumption. It postulates that, when the locally 
observed flow coincides with the one associated to the aggregate path preference tt, 
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drivers choose to join link e € with frequency equal to the ratio between the flow 
fl and the total outgoing flow ^2j e£ + fj- 

Assumption 6 instead models the drivers' myopic behavior in response to varia- 
tions of the local congestion levels. It postulates that, if the congestion on one link 
increases while the congestion on the other links outgoing from the same node is kept 
constant, the frequency with which each of the other outgoing links is chosen does not 
decrease. It is worth observing that Assumption 6 is reminiscent of Hirsch's notion of 
cooperative dynamical system [9, 10]. 

Example 3. An example of local decision function G v satisfying Assumptions 5 
and 6 is the i-logit function. The i-logit route choice with sensitivity 7 > is given by 

r v (f x fe exp(-7(/ e - fj)) 

Ge {f£ * ' n) = E^ + /;cxp(- 7( /,-/;)) • ^ 

for every e G £+, < v < n. 

For every non-destination node v £ {0, 1, . . . , n — 1}, and outgoing link e e £+, 
conservation of mass implies that 

d _ / Gl(f £ +,n)-f e if v = 



dt p e -H e (f,Tr), H °V> n >-\ (£. f . e -f j )Gv(fs+>*)-fe if l<v<n. 

(2.10) 

2.3. Objective of the paper and main result. The objective of this paper 
is to study the evolution of the coupled dynamics 



^ = v (F h (f)-n) 



(2.11) 



where F h is the perturbed best response function defined in (2.4), 77 > is the rate 
at which global information becomes available, H(f,n) = {H e (f,n) : e e £}, with 
H e defined in (2.10), and / and p are related by the functional dependence (2.1). In 
particular, our analysis will focus on the double limiting case of small 77 and small h. 
We shall prove that, in such limiting regime, the long-time behavior of the system is 
approximately at Wardrop equilibrium [20, 15]. The latter is a configuration in which 
the delay is the same on all the paths chosen by a nonzero fraction of the drivers. 
More formally, one has the following: 

Definition 2.2 (Wardrop Equilibrium). An admissible flow vector f w e T is 
a Wardrop equilibrium if f w = An for some w € II such that, for all p G V , 

tt p >0 => {A'T(An)) p < (A'T(Aw)) q , VqeV. (2.12) 



Existence and uniqueness of a Wardrop equilibrium are guaranteed by the following 
standard result: 

Proposition 2.3 (Existence and uniqueness of Wardrop equilibrium). Let As- 
sumptions 1-3 be satisfied. Then, there exists a unique Wardrop equilibrium f w e 
T. 

Proof. It follows from Assumption 2 that, for every e G £, the delay function 
T e (/ e ) is continuous, strictly increasing, and such that T e (0) > 0. The proposition 
then follows by applying Theorems 2.4 and 2.5 from [15]. ■ 
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The following is the main result of this paper. It will be proved in Sect. 3 using 
a singular perturbation approach. 

Theorem 2.4. Let Assumptions 1-6 be satisfied. Then, for every initial condition 
7r(0) £ int(<S('P)), p(0) £ (0,+oo) £ , there exists a unique solution of (2.11). Moreover, 
there exists a perturbed equilibrium flow f h £ J- such that, for all rj > 0, 

limsup||/(i)-/ /l || <S( V ), (2.13) 

t— ^ + 00 

where S(i]) is a nonnegative-real-valued, nondecreasing function of n > 0, such that 
lim^o S(rj) = 0. Moreover, for every sequence of admissible perturbations {hk ■ k £ 
N} such that limfc( +oc ||/ife||oo = 0, and linifc^ +00 Uh k = II, one has 

lim f hk = f w . (2.14) 



Theorem 2.4 states that, in large time limit, the flow vector /(f) approaches a 
neighborhood of the Wardrop equilibrium, whose size vanishes as both the time-scale 
ratio r] and the perturbation norm | \ h\ \ QO vanish. While a qualitatively similar result is 
known to hold [18] in a classic evolutionary game theoretic framework (i.e., neglecting 
the traffic dynamics, and assuming it is instantaneously equilibrated, as in the ODE 
system (2.7)), the significance of the above is to show that an approximate Wardrop 
equilibrium configuration is expected to emerge also in our more realistic model of two- 
time scale dynamics. Therefore, our results provide a stronger evidence in support of 
the significance of Wardrop's postulate of equilibrium for a transportation network. 
In fact, they may be read as a sort of robustness of such equilibrium notion with 
respect to non-persistent perturbations. 

3. Proofs. In this section, Theorem 2.4 is proved. First, observe that, thanks to 
the continuous differentiability of F h , G v , and p, standard analytical arguments imply 
the existence and uniqueness of a solution of the the initial value problem associated 
to the system (2.11), with initial condition p(0) £ (0, +oo) £ , 7r(0) e int(n). 

In order to prove the rest of the statement, we shall adopt a singular perturba- 
tion approach (e.g., see [13]), viewing the traffic density p (or, equivalently, the traffic 
flow /) as a fast transient, and the aggregate path preference ir as a slow component. 
Hence, we shall first think of 7r as quasi-static (i.e., 'almost a constant') while analyz- 
ing the fast-scale dynamics (2.10), and then assume that / is 'almost equilibrated', 
i.e. close to f n , and study the slow-scale dynamics (2.5) as a perturbation of (2.7). 
We shall proceed by proving a series of intermediate technical results, gathered in the 
following subsections. 

Before proceeding, we introduce some notation to be used throughout the section. 

Let 

Pe : = t i e 1 (.fe) , *e 1= Sgn (p e - p%) = Sgn (f e - ft) 

denote, respectively, the density corresponding to the flow associated to the path 
preference n, and the sign of the difference between it and the actual density p e . 
Finally, fix some a € (0, 1), and define 

n— 1 n— 1 

V(f,n):=J2* V E l/e-/ri. W(p,n):=J2 aV E \Pe ~ PlV 

v =° ee£+ v=0 ee£ + 
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3.1. Stability of the fast-scale dynamics. We gather here a few properties of 
the fast-scale dynamics. Our results will essentially amount to showing that V( f, n) 
and W(p, 7r) are Lyapunov functions for the fast-scale dynamics (2.10) with stationary 
path preference n. 

The following result is a consequence of Assumptions 5 and 6 on the drivers' local 
decision function. 

Lemma 3.1. For all tt E IT, v E {0, . . . , n - I}, and f £ + E T v , 

]T o e (KGlU £t ,v)-tt) <o, 

where A" := £ ee£ + fe- 

Proof. Throughout this proof, the explicit dependence of G v e on 7r will be dropped. 
Define J := {e E £+ : f e > K := {e E £+ : / e < /J}, and let Gj := E j£l7 G), 
G K. ■= J2keic G fc> and G J C '■= J2ee£+\J G e- First > observe that, since £ ee£ + G e = 1, 
one has that VGj = -VGjc. Now, we are going to show that 

Gj(%+)-Gj(f £ +)>0, (3.1) 

by writing the difference above as a path integral of WGj( ■ ) first along the segment 

£ + 

Sj from f £ + to the point /* E R + " with f* := f?, for j E J and /* := f e for 
e e £+ \ and then along the segment S)e from /* to /*■. In this way, one gets: 



Gj{fl t )-Gj{f £t ) = f VGj(f £ +)-df £ ++ f WGjCfei)^hi 

J Sj J Sic 

= - ! VG JC (4 + ) • df £ + + f VGj (f £ + ) • d/ £ + . 

J S.t J Sir. 



(3.2) 

Assumption 6 implies that dGj^/dpj > for all j E J, and dGj/dpu > for all 
k E JC. In turn, this implies that VGjc • df £ + < along Sj, and VGj • df £ + > 
along Sjt- This and (3.2) prove (3.1). In a very similar fashion, one proves that 

GK{f £ +)-G K {r £ +)>Q. (3.3) 

Now, observe that Assumption 5 implies that AJGg(/^+, n) — f£. From this, (3.1), 
and (3.3), it follows that 

> K(Gj{f £ +)-Gj{fl t ))-K(G K (f £ +)-G K {r £t )) 

= X! CTe (K G e(f£+) ~ KGl(f £ +)j 
= ^ G e (K G e(fe+) - fe) . 

which proves the claim. ■ 

We now proceed to analyzing, for a fixed global decision it £ IT, the fast scale 
dynamics (2.10). Let 

V V + (L n) := Y, ee£ + \fe ~fe\ , v = 0, 1, . . . , n - 1 , 
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be the ^-distance between the current flows on the outgoing links of v, and the flow 
associated to the aggregate path preference ir, and 



V-(f,7r):=\X:-\-\ , « = l,2,...,n, 

with AJ := J2 e e£ + fe anc ^ K> '■= J2 e e£~ ^ e tne absolute difference between the 
current flow incoming in node v, and the one associated to the aggregate path pref- 
erence 7r. Also, let V^~(f, 7r) := 0. 

Lemma 3.2. For all v = 0, 1, . . . , n - 1, n e IT, and f e J", 

E CT e^e(/,7r) < -V+{f,n) + V-(f,n) . 

ee£ + 

Proof. Writing G v e for G^(f £ + 7 n), and using Lemma 3.1, one gets that 
£a e #e(/,7r) = E + ^(A-G^-/ e ) 

= 6 E ^e(A" - A-)G- + £ <7 e - /T) + £ <7 e (/J - /e) 

< ^ - K\ - E e£ 5+ I/? - /el 

= -K + G/» + Kr(J», 

which proves the claim. ■ 

By combining Lemma 3.2, and Assumption 1, one gets the result below. Recall 
that we are using the convention d|a;|/da; = sgn(a;), for all x € R. 
Lemma 3.3. For every f = p(p) e J ' , and n e II, 

V p W(p, n)'H(f, tt) < -(1 - a)V(f, n) . 



Proof. Observe that, thanks to the acyclicity of the graph as per Assumption 1, 
if e e £~ n for some nodes v and w, then necessarily v > w + 1. Since a < 1, it 
follows that 

a%-(e)%( e )<^+%-(e)i £ +(e), 

for every 1 < u < n, and < w < n — 1. Hence, 

E ^k-c/» < E E a "i/e-/n 

0<v<n 0<v<7i e ££~ 

= E E E at, v( e )%( e )i/ e -^i 

1 < v < ?i 0<it)<n e££ 

< E ^E^^i^-/:! E W e ) 

< id < n e££ l<i)<n 

< a E E l^-^l 

0<w<ri e££ + 

= aV(f,n), 
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where the last inequality follows from the fact that J2v=i ^-£~ ( e ) — S"=i ^st ( e ) = 1- 
Thus, Lemma 3.2 implies that 



V p W(p,7r)'H(f,n) = £ a v ]T o- e H e (f,n) 

0<v<n ee £+ 

0<v<n 0<v<n 

< aV(f,w)-V{f,ir), 



which proves the claim. 



3.2. Boundedness of the traffic densities. We shall now prove a couple of 
results guaranteeing that the traffic density on every link e e £ remains bounded in 
time. We start with the following result, guaranteeing that, on every link e € £, the 
flow associated to the current path preference, f£(t), stays eventually bounded away 
from the maximum flow capacity C e . Its proof relies on Assumption 4. Recall that 
our formulation allows for both the cases of finite and infinite maximum flow capacity 
on a link. 

Lemma 3.4. For every admissible perturbation h, there exists t € M+ ; and, for 
every link e G £, a positive finite constant C e , dependent on h but not on n, such 
that, for every initial condition 7r(0) € int (S(P)), p(0) G (0, +oo) £ , 

fZ(t)<C e <C e , Vt>t , Vee£. 

Proof. The fact that f£(t) < 1 for all e <G £ follows from the fact that the arrival 
rate at the origin is unitary. Therefore, for all e e £ with C e > 1 or C e — oo in 
particular, the lemma follows trivially with C e = 1 and to = 0. We now prove the 
lemma for all e £ £ with C e < oo. Recall that, by Assumption 4, the domain of the 
admissible perturbation h is a closed set 11^ C int (IT). This, in particular implies 
that, for all e e £ with C e < oo, 

n e := C e - sup{(yLr) e : tt e U h } > . 

For every link e € £ with C e < oo, it follows from (2.4) that 

C e -n e = sup{(A-7f) e : 7r e U h } 

> sup{(Aargmin{7f' J 4'T(/) + h(n) : n e II fc }) e :/eJ} (3.4) 
= sup{(A^(/)) e : feF}. 

Hence, one has for every link e £ £ with C e < oo, 

A#( t ) = „ (A(F"(/(t)) - 7r(t))) e < ry (C e - Ke - . 
Then, Gronwall's inequality implies that 

/T(*) < /;(0)e- r " + (1 - e-"*)(C e - Ke ) < e"* + C e - Ke , 

for alH > and for every link e G £ with C e < oo. The lemma for e € £ with C e < oo 
now follows from the above, by choosing, e.g., C e := C e — k with « := ^ min{« e : e e 
£ s.t. C e < oo}, and i : = — V 1 logK. I 
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The following result shows that the actual flow f e (t) also stays bounded away 
from the maximum flow capacity C e . 

Lemma 3.5. For every admissible perturbation h, there exists positive finite con- 
stants rj* , and C e , for every e <G £, such that, for every rj < n* , and every initial 
condition tt(0) G int(S(7>)) } p(0) G (0,+oo) £ , 

f e (t)<C e <C e , Vt>0, VeGf. 



Proof. Let ((t) := W(p(t),n(t)), and x(t) ■= V(f(t),n(t)). Observe that, thanks 
to Lemma 3.4, there exists a positive real constant C e for every e G £, and to > 0, 
such that, for every t>t , 

Pl{t)<p* e , p* e :=ii-\C e ), Vee£. (3.5) 

Since p1(t) > 0, the above implies that, if \p e (t) — P^it)\ > 2p* for some e G £ 
and t > t , then necessarily for that e G £, we have that f e (t) — f£(t) > xt, where 
X* := p e (2p*e)-C e > 0. Now, let (* := 2\£\ max{p* : e G £}, and X * ■= a™ -1 min{x* : 
e e £}. Note that W{p,n) < \£\ max{|p e - p£)| : e G £}, and V{f,n) > a™ -1 ^ - /"I 
for every e G £. Therefore, it follows that, for any t > t , if ((t) > (* , then for 
some e' G £, we have that \p e '(t) — pj| — ^Pe' f° r * — *o- This in turn implies that 
X(t) > Xe' — X*- Therefore, in summary 

C(*)>C* => xW>x*, vt>t . (3.6) 

On the other hand, observe that (3.5) implies that there exists some £ > such 

that 

1 7 » < ^ , Vt > t • 
By combining the above with Lemma 3.3, one finds that, for any u,t> to, 

< / V p W(p, 7 r) / J ff(/, 7 r)ds 

A r t (3.7) 
+ / E«"E T^T \iAF\r))e (An) e \ds 

Ju 0<v<n ee£ + VekPe) 



< 



f (-(l-a)x(s)+2 V l)ds. 

J u 



Now, define rj* := (1 — a)x* /(%£)■ Assume, by contradiction, that limsup^^ 
fe(t) > Ce for some e € £■ Since / e (i) = p e (p e (t)) < C e for every i > 0, this implies 
that lim sup t _j. +00 p e (t) = +oo. This, together with (3.5) implies that limsup t ^ +co 
((t) — +oo. Then, in particular, the set T := {t > : ((i) > C( s ) jVs < £} is an 
unbounded union of open intervals, with linii^7 _ ^_j._|_ 00 

C(t) = +oo. This, and (3.6) 
imply that there exists a non negative constant t* > to such that 



x(t) >x% Vt > t* . 
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For every r\ < 77*, Equation (3.7) and the above give 



C(t)-C(«)= / (-(l-a)x(s)+2 V l)ds< I (-(l-a)x*+2^)ds<0 

J u Ju 

for every t > u > t* such that i and u belong to the same connected component 
of T. But this contradicts the definition of the set T. Hence, if 77 < 77* , then 
lim sup t _j. +00 f e (t) < C e for every e € £. Since on every compact time interval I C 
[0, +00), one has sup t6l / e (t) = f e (i) < C e for some i G I, the foregoing implies the 
lemma. ■ 



The result below is a consequence of Lemma 3.5, and will prove useful in the 
sequel. 

Proposition 3.6. There exists K > ; and t\ > such that, for every initial 
condition tt(0) e int {S{V)), p(0) G (0, +oo) f , HV^tt^))!! < K for all t > h. 

Proof. First, observe that, thanks to Lemma 3.5, there exists T* > such that 
< T* for all t > 0. Thanks to this, and Assumption 4, one has that 
F (f(t)) G int(n h ), and V vh{F h {f (t))) = -<f>A'T(f(t)), where recall that $ = 
J — I? 13 1" 1 11' is the projection matrix corresponding to the projected gradient with 
respect to tt on S(P). Hence, \\V^h(F h (f(t)))\\ < ||$||P'||T*, which implies that 
there exists a convex compact JC C int(n;j) such that F h (f(tj) e JC for all t > 0. 
Define 

/•t 

A(i) := 77 (l - e^') -1 / e-^ t -^F h (f(s))ds. 
Jo 

As A(i) is an average of elements of the convex set JC, necessarily A(t) e JC for all 
t > 0. Then, tt(<) = e'^^O) + (1 - e _,7 *)A(/;) approaches /C, which implies that, 
for large enough t, ir(t) £ K\ C int(n^), where /Ci is a closed subset of int(n/j) that 
contains JC. Hence, after large enough t, say t\, V„-/i(7r(/j)) stays bounded. ■ 



3.3. Estimating the distance between the current density and the one 
associated to the current path preference. We analyze here the behavior in time 
of W(p(t),ir(t)). First, we have the following result, characterizing the variation of 
W(p, 7r) as a function of tt. Recall that = (I — |7 , | _1 ll')V 7r denotes the projected 
gradient with respect to n on S(V). 

Lemma 3.7. There exists I > 0, and t > 0, such that, for every initial condition 
tt(0) e mt(S(V)), p(0) € (0,+oo) £ , 

0/ 

v v w(p(t)Mt))'(F h m) - At)) < . *o ■ 



Proof. First, observe that, thanks to Lemma 3.4, one has that there exists to > 
such that l e := sup{l / p' e (p^ (t)) : t > t } < +00. Put I := m&x{l e : e G £}. Then, 
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for every path p € V, and every t > to, one has 



dW{p,ir) 



E «' E * 



0<t)<n 



E a " E ^^^'(EpVp) 

1 



0<t)<n 



< E «' E A * 



< 



< 



0<v<n ee£ + ^e\l J e) 
0<v<n ee g + 



1 



(3.8) 



where the third inequality follows from the fact that, thanks to Assumption 1 on the 
acyclicity of the network, each path peP passes through at most one link e e £+. 
Therefore, 



21 

1 - a 



> 



d-K, 



W(p,7T) 



+E^ 



J)_ 



W(p,n) 



p y p y 

= V w W(p,Tr)'(F h (f)-iT), 



where the first inequality follows upon recalling that both F h {f), and it are probability 
vectors over the path set V, and by using (3.8). ■ 

We can now combine Lemmas 3.3 and 3.7, in order to get the following estimate 
of the behavior in time of W(p(t),ir(t)). 

Lemma 3.8. There exist I > 0, L > 0, n* > and t > such that, for every 
initial condition n(0) e int(S(V)), p(0) € (0, +oo) £ , 



2r]lL 



< 



(I -af 



W(p(t ),7T(t )) - 



2nlL 

(T^p 



e -V(*-*o) y t >t Q , V77 > 77*. 



Proof. Define x(t) := W(p(t),n(t)). Notice that, thanks to Lemmas 3.4 and 
3.5, there exist L > 0, i]* > and to > 0, such that, for any 77 < 77*, \p e (t) — 
PeWI — L\fe{t) — fe(t)\ f° r every e e £, t > to- This in particular implies that 
V(f (t),n(t)) > W(p(t),ir(t))/L = x(t)/L, for all rj < rj* and t > t . Observe that 
W(p, 7r) is a Lipschitz function of p and ir, while both p{t) and Tr(t) are Lipschitz 
on every compact time interval. Therefore, x(t) is Lipschitz on every compact time 
interval, and thus differentiable for almost every t > to- For every t at which x(t) is 
diffcrcntiable, Lemmas 3.3 and 3.7 imply that 

d ;x(t) < V p W(p,n)'H(f,7r) + r 1 V 7r W(p,iry(F h (f)-n) 

2nl 



dt 



< -(l- a )V(f,ir) , . 
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< 



Then, the claim follows from Gronwall's inequality. ■ 

3.4. Proof of Theorem 2.4. We now proceed to proving Theorem 2.4. Con- 
sider the following candidate Lyapunov function: 

ffe 

^(tt) := V / T e (s) ds + h(n) . (3.9) 
eee Jo 

Since T e (/ e ) is increasing, one has that each term J Q e T e (f e )df e is convex in /J. 
Hence, the composition with the linear map tt \-t f* = J2 p ev ^-&p^p is convex in tt. 
Since h(ir) is strictly convex by assumption, one gets that ip h (ir) is strictly convex as 
well. Therefore, tp (tt) admits a unique minimizer 

TT h := argmin {tp h {n) : tt <G H h } . (3.10) 

Let f h := Aw h . Then, we have the following: 

Lemma 3.9. Let {h k ■ k <G N} be any sequence of admissible perturbation func- 
tions such that lim ||ftfe||oo = 0; limlL^ = II. Then, 

k k 

lim f h » = f w . 

k— >-\-oo 

Proof. Write n k for ir h « , F k for F h " , and n fc for . Since {Air k } C AU 7 and AU 
is compact, there exists a converging sub-sequence {Air ki : j g N}. Let us denote by 
/* := limj Aiz k ' g All its limit, and choose some tt* g U such that /* = Air*. Notice 
that, since sup{T e (/ ( I r ) : tt £ H^} < +oo, Assumption 4 implies that the minimizer 
in (3.10) has to be in the interior of 11^. As a consequence, one finds that necessarily 
V n h(w k i) = -<&A'T(ATT k i), which in turn implies that F k ? (Air k i) = TT k ' . Then, using 
(2.4), one finds that 

{Air k > )'T(Air k i ) + h kj (tt^ ) < {Att)'T{Att^ ) + h kj (tt) , (3.11) 

k — 

for all tt g IL^. . Now, fix any tt G II. Since life — > II, one has that tt € 11^ for all 
sufficiently large values of j. Hence, passing to the limit as j — > +oo in (3.11), one 
finds that 

(tt* )'A'T{Att* ) < tt'A'T{Att* ) , Vtt g n . 

In turn, the above can be easily shown to be equivalent to the condition (2.12) char- 
acterizing Wardrop equilibria. From the uniqueness of the Wardrop equilibrium, it 
follows that necessarily /* = f w . Then the claim follows from the arbitrariness of 
the accumulation point /* , and the compactness of An. ■ 

We shall now estimate the time derivative of ip h {iT) along trajectories of our 
dynamical system. For this, define 

r(t) :=^(7r(t)). 

Then, one has 

= V a'(F h (f(t))-TT(t)) (3-12) 

= na' (F h (An(t)) tt(*)) + r,a' (F h (f(t)) F h (ATr(t))) . 
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where a := &A'T(An(t)) + V n h(ir{t)). 

Lemma 3.8 implies that there exists t 2 > 0, 77* > and Mj > such that, for any 
r] < rj*, W(p(t),ir(t)) < r\M\ for all t > t 2 . From the definition of W, it also follows 
that W(p, 7r) > p n \\i for all p,7r. Moreover, following Assumption 2, with 

L := max e£ £ §^(0), we also have that ||/ — Att\\i < L||p — p^Hi for all / = p(p) and 
7r. Combining all these relationships, one can see that there exists a M > such that, 
for any i] < r/* , 

\\f(t)-Air(t)\\<r)M, Vt > t 2 . (3.13) 

Moreover, recall that F h is diffcrentiable on J 7 , and that, thanks to Lemmas 3.4 and 
3.5, for rj < i]* , both f(t) and Aw(t) are eventually confined in a compact JC C J 7 . 
This implies that 

||F"(/(t))-F fc (A7r(t))||<^ 

for some positive constant ^i, 77 < ry* and sufficiently large values of t. On the 
other hand, Lemma 3.4 and Proposition 3.6 imply that both T(An(t)) and V ff /i(7r(t)) 
are eventually bounded, so that ||a|| < K 2 , for some positive constant K 2 and large 
enough t. It follows that the second addend in the last line of (3.12) can be bounded 
as 

V a' (F h (f(t))-F h (An(t))) < K V 2 , Vt? < rf , Vt>t 3 , (3.14) 

for some sufficiently large but finite value of £3, where K — K\K 2 . Now, observe that, 
for every ir, 

<$>A'T{Air) = -V 7I h(F h (An)) , 
so that the first addend in the last line of (3.12) may be rewritten as 

a'(F h (An(t))-n(t))=-r(n(t)), (3.15) 

where 

T(tt) := (v„h{F h {Ait)) - V w /i(0)' (F h (An) - tt) . 
It follows from (3.12), (3.14), and (3.15), that, for 77 < rj* and t > t 3 , 

fr(t) < -vr(n(t)) + M V 2 . (3.16) 

From the strict convexity of h(n) on the simplex II, one finds that T(ir) > for all 
tt, with equality iff n = ir h . Now, let 

( su P {\\Att - f h \\: T(ir) < Mx} + Mx if < x < r/* 
[X > : ~ \ Cj\t\ if x > rj* , 

where C := max{C* e : e e £ }, with C e as defined in Lemma 3.5. It can be verified 
that S(x) is right-continuous, nondecreasing, and such that 3(0) = 0. Then, (3.13) 
and (3.16) imply that, for r\ < rf , 

limsup||/(t)-/' l || <S( V ). 

For T] > rj* , the above is clearly true since f(t) € [0,(7] by Lemma 3.5 and f h E 
AH C [0, 1] £ . Together with Lemma 3.9, this completes the proof of Theorem 2.4. 
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Fig. 4.1. The graph topology used in simulations. 



4. Simulations. In this section, we present results from numerical experiments. 
We performed several experiments with different graph topologies and for values of 
r\ ranging from 0.01 to 100. In all the cases, we found that the trajectories con- 
verge exactly to the perturbed Wardrop equilibrium, i.e., 8(r\) in Theorem 2.4 was 
estimated to be uniformly zero. We suspect that this might be because of the expo- 
nential convergence also of the slow scale dynamics. Additionally, we compared the 
convergence of the trajectories corresponding to local decision function from Example 
3 with trajectories corresponding to local decision function of the form 

G l (fet > *) = fel E fi > V 4+ e Ve e # ■ ( 41 ) 

The latter corresponds to the case when the drivers do not take into account the local 
observation on the currently observed flow, and always act in a way that is consistent 
with their aggregate path preference. We found that the trajectories corresponding 
to local decision function in (4.1) converged faster than the trajectories corresponding 
to the local decision function in Example 3. 

We demonstrate these findings through an illustrative example. For this example, 
the parameters were selected as follows: 

• graph topology Q as shown in Figure 4.1, 

• link- wise flow functions as given by (2.3) with C\ — 2 and 9 e = 1, for all 

e € £; 

• F h as in (2.6) with /3 = 1, 

• G as in (2.9) with 7 = 1, 

• initial conditions: 7r e (0) = 1/15 for all e G £, p ei (0) — p £l2 (0) = 5, p e2 (0) = 
Pe6 (0) = p fi8 (0) = 7, Pe M = p e7 (0) = 3, p e4 (0) = 6, p e , (0) = 1, p eg (0) = 9, 
PeiB (0) = 10, p ei3 (0) = 12, Pei4 (0) = 4, p eiB (0) = 8. 

• 7/ = 0.1. 

For these values, p h :— p~ 1 {f h ) was numerically calculated by implementing a gradi- 
ent descent algorithm for the potential function as given in (3.9). The evolution of the 
1-norm distance of p from p h is plotted on a log-linear scale in Figure 4.2 for two cases: 
(i) local route choice decision function of Example 3, and (ii) local decision function 
given in (4.1). Figure 4.2 also shows that there is no significant difference between 
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Fig. 4.2. Log-linear plot for comparison of the evolution of \\p{t) — p h \\i for the local decision 
function of Example 3 versus the local decision function of (4-1). 



the convergence of trajectory corresponding to local decision function in (4.1) and the 
trajectory corresponding to the local decision function of Example 3. However, as we 
increase rj, we observed that the trajectory corresponding to the local decision func- 
tion in (4.1) converge faster than the trajectory corresponding to the local decision 
function of Example 3. 

5. Conclusion. In this paper, we analyzed the stability of Wardrop equilibria 
in dynamical transportation networks characterized by dual temporal and spatial 
scales of the drivers' route choice behavior. We showed that, if the frequency of 
updates of path preferences is sufficiently small, then the state of the transportation 
network ultimately approaches a neighborhood of the Wardrop equilibrium. The 
technical approach relied on establishing relevant properties for the resultant two 
time-scale dynamics independently using tools from evolutionary game dynamics, and 
cooperative dynamical systems, and then using singular perturbation techniques to 
establish sufficient conditions for the stability of the Wardrop equilibrium for the 
coupled system. Our results contribute to providing a stronger evidence in support of 
the significance of Wardrop's postulate of equilibrium for a transportation network. 
They may be read as a sort of robustness of such equilibrium notion with respect to 
non-persistent perturbations of the network. 

There are several possible directions for future work. We plan to formally justify 
our dynamical model as a macroscopic approximation of the underlying driver level 
microscopic process. We also plan to extend our analysis to the case with multiple 
origin-destination pairs and possibly cyclic topologies. We also plan to study the 
effect of persistent, and possibly adversarial, perturbations on the traffic dynamics 
under driver behavior model similar to the one considered in this paper, e.g., see [5]. 
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